Doping dependence of the superconducting gap in Tl 2 Ba 2 Cu06+5 from heat transport 
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We present low-temperature thermal conductivity measurements on the cuprate Tl2Ba2Cu06+5 
throughout the overdoped regime. In the T — > limit, the thermal conductivity due to d-wave 
nodal quasiparticles provides a bulk measurement of the superconducting gap, A. We find A to 
decrease with increasing doping, with a magnitude consistent with spectroscopic measurements 
(photoemission and tunneling). This argues for a pure and simple d-wave superconducting state in 
the overdoped region of the phase diagram, which appears to extend into the underdoped regime 
down to a hole concentration of ~ 0.1 hole/Cu. As hole concentration is decreased, the gap-to- 
T c ratio increases, showing that the suppression of the superconducting transition temperature T c 
(relative to the gap) begins in the overdoped regime. 

PACS numbers: 74.25.Fy, 74.72. Jt 
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Over the past decade, low-temperature thermal con- 
ductivity has emerged as a powerful probe of the low- 
energy quasiparticle excitations in the cuprates. Unlike 
conventional s-wave superconductors, the nodes in the d- 
wave superconducting order parameter of cuprates lead 
to a finite quasiparticle density of states at zero energy re- 
sulting in "universal" transport in the low-temperature 
limit (T — > 0). 1,2 By performing measurements in the 
T — ► limit, the complicated physics of quasiparticle 
transport in these materials is greatly simplified and 
meaningful quantitative comparisons between mean-field 
theories in the self-consistent T-matrix approximation 
(SCTMA)' 3 .!-^ and experiments have been made^i^ 

In the mean-field theory, the low-T thermal conduc- 
tivity is proportional to the slope of the gap near the 
nodes, SA/S(j)\^ w u. By assuming a gap with a simple 
d-wave form (A = Ao cos(2</>)), the gap maximum, Ao, 
can be inferred from SA/S<p, making thermal conductiv- 
ity a measure of the superconducting gap. This measure 
is notable for a number of reasons that make it an ex- 
cellent complement to more conventional, spectroscopic 
measurements of the gap. First, it is a very low energy 
measurement of the gap, the energy scale being set by 
temperature (and impurity bandwidth, 7). Second, it 
measures the gap near the nodes (in the T — > (it, w) 
direction) as opposed to the anti-nodes (irfi). Finally, 
thermal conductivity provides a bulk measure of the gap, 
unlike surface sensitive techniques such as photoemission 
(ARPES) and tunneling. Together these properties al- 
low one to impose important constraints on the nature 
of superconductivity in cuprates. 

The cuprate Tl 2 Ba 2 Cu0 6 +5 (Tl-2201) is an ideal ma- 
terial for quantitative investigations of d-wave super- 
conductivity in the cuprates. It is a single-layer mate- 



rial with T™ ax ~ 90 K, comparable to the well-studied 
cuprates Bi 2 Ba 2 CaCu 2 08+,5 (Bi2212) and YBa 2 Cu30 y 
(YBCO). By varying the oxygen concentration, the ma- 
terial can be doped from optimal doping far into the over- 
doped region of the phase diagram. In this region, the 
metallic state is well-characterised in terms of a single 
coherent Fermi surface and there are no indications of 
competing phases, such as encountered in La^zSr^CuO^ 
(LSCO), for example. Moreover, Tl-2201 can be prepared 
with reasonably low levels of disorder because the dopant 
atoms are interstitial oxygen residing well away from the 
Cu0 2 planes, instead of the more obtrusive cation dop- 
ing located close to the Cu0 2 planes (as in LSCO, for 
example) . 

In this paper, we present low-temperature thermal con- 
ductivity measurements of Tl-2201 as a function of dop- 
ing, spanning the overdoped region of the phase diagram. 
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Expanding upon previous work on cuprates 
analyze our results using the mean-field SCTMA theory 
to provide a determination of the superconducting gap 
from thermal conductivity. We find agreement between 
thermal conductivity and spectroscopic measurements of 
the gap, confirming both the d-wave form of the super- 
conducting gap as well as the validity of the mean-field 
description of low-temperature thermal conductivity in 
overdoped T12201. We track the gap as a function of 
doping throughout the overdoped region and find a gap- 
to- T c ratio, /1 = Ao/fcsT c , that increases with decreasing 
doping, starting with a value close to the weak-coupling 
BCS value (/.i = 2.14) at high doping and reaching a 
value twice as large near optimal doping (/1 = 5 — 6). 
This shows that the suppression of the onset of coherent 
superconductivity (or T c ) so characteristic of the under- 
doped regime actually begins above optimal doping. 
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FIG. 1: (colour online) Thermal conductivity of one Tl-2201 
sample, measured at two different hole concentrations, as in- 
dicated by the respective T c values, plotted as k/T vs T below 
1 K. The lines are fits to Eq. [2] A three-fold decrease in T c 
corresponds to a six-fold increase in the residual linear term 
(which is inversely proportional to the superconducting gap, 
via Eq. [5}, hence a two-fold decrease in the gap-to- T c ratio 
with doping. 



Five flux-grown single crystals of T^E^CuOg+s are 
used in this study£ii£ The platelet samples have approx- 
imate length and width of ~ 0.5 mm in the ab plane and 
thickness x 20-65 jxm along the c axis. The heat current 
is applied in the ab plane and the magnetic field along 
the c-axis. The concentration of holes in the CuC>2 plane, 
p, is set by tuning the excess oxygen content in the TIO2 
layer. Concentrations between optimal doping and heav- 
ily overdoped were achieved by annealing the samples in 
atmospheres ranging from vacuum (~5x 10~ 7 torr) at 
500 °C near optimal doping to 50-60 mbar oxygen at 
400-480 °C for heavily overdoped samples. The value of 
p for each sample was estimated from the measured T c 
(defined as the temperature at which the electrical resis- 
tivity has fallen to zero), using the well-known empirical 
formulaii 



l-82.6(p-0.16) 2 



(1) 



where T" lax = 90 K is taken as the nominal T c of 
optimally-doped Tl-2201. The value of T c and p for each 
sample is listed in Table I. 

Thermal conductivity, n, was measured down to 50 mK 
in a dilution refrigerator using a standard 4- wire steady- 
state technique described elsewhere.— Electrical contacts 
to the samples were made using evaporated gold pads, 
which were annealed to improve the contact conductance 
in the same anneal used to set the oxygen content. Silver 
wires were then attached to the contact pads using silver 
paint. 
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FIG. 2: (colour online) Thermal conductivity of four Tl-2201 
samples, each with a different hole concentration, as indicated 
by their respective T c values, plotted as k/T vs T below 1 K. 
The lines are fits to Eq. [2] 



TABLE I: Residual linear term, Ko/T, in the thermal conduc- 
tivity of Tl-2201, from fits of k/T vs T to Eq.EJ for five single 
crystals with hole concentration p estimated from measured 
T c . The gap maximum, Ao, is evaluated from ko/T using 
Eq.[5] The error applies to Ko/T, Aq and Ao/fc_eT c . 
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"Denotes the same sample measured at two dopings. 



The largest error in the magnitude of k, ranging from 
15 to 25 % (see Table I), comes from the measurement of 
sample dimensions for the appropriate geometric factor. 
In order to eliminate this dominant source of uncertainty 
when comparing two different hole concentrations, one 
of the samples was measured twice with the same con- 
tacts, at two different dopings, T c = 84 K and 26 K (see 
data in Fig. 1). This was achieved by first evaporating 
gold pads onto the sample and annealing the sample to 
set the doping, as described above. After measuring the 
sample, the wires and silver paint were removed, the sam- 
ple was reannealed to a different doping and subsequently 
re-measured with the gold pads still in place. 

The thermal conductivity measured below 1 K is shown 
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in Fig. [T] for this sample and in Fig. [2] for four other sam- 
ples, each at a different doping. To extract the electronic 
contribution to k the data are fit to 



- = f + AT, 



(2) 



where kq/T (intercept) is the electronic contribution 
to the total thermal conductivity and the second term 
(slope) is the phonon contribution, n p h/T = AT. The 
values of kq /T from fits to Eq. [5] are listed in Table I, 
along with the error bars for each sample, estimated from 
the statistical error in the fits (< 1.5%) and the error in 
the geometric factor (15 — 25%). These results can be 
compared to a previous study on a strongly overdoped 
Tl-2201 crystal^ with T c = 15 K, which gave k q /T = 1.4 
mW K~ 2 cm -1 and an excellent fit to Eq. O 

Note that in these overdoped samples n p h oc T 2 (over 
a decade in temperature). This is in contrast to previ- 
ous reports of the low-temperature thermal conductivity 
of cuprates, measured typically at or below optimal dop- 
ing, where n p h cx T 25 ~ 3 , consistent with phonons being 
scattered from the boundaries of the sample^^^ ' 12 ^ 3 
In overdoped Tl-2201, with its higher carrier concentra- 
tion, we attribute the T 2 dependence of K p h to dominant 
electron-phonon scattering, known to yield K p h cx T 2 at 
low T in a metal^ Smith has shown that this is indeed 
the theoretically expected dependence of phonon trans- 
port in a d-wave superconductor^ in the range 0.1 - 
1.0 K. We emphasize that the main results of this pa- 
per, which stem from the doping dependence of kq/T, 
are unaffected by alternative methods to fit K p h, such as 
a power-law with adjustable power (k/T = a + bT 1 ). 

The T-linear electronic contribution to k is due to 
transport from d-wave nodal quasiparticles in the limit of 
ksT -C7< Ao, where 7 is the impurity band width.— 
In this limit, the residual linear term in the thermal con- 
ductivity is universal, in the sense that it is independent 
of scattering rate, and given by^ 
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where n is the number of CuC>2 planes per unit cell 
(= 2 for Tl-2201) and c is the c-axis lattice constant 
(= 23.2 A for Tl-2201). The ratio v F /v A is the ratio 
of quasiparticle velocities normal and tangential to the 
Fermi surface at the node, respectively given by the Fermi 
velocity in the (0,0) - (-7T, 7r) direction, vf , and 
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FIG. 3: (colour online) Upper panel: The gap maximum, Ao, 
vs hole concentration, p, derived from thermal conductivity 
measurements on Tl-2201 (circles) and YBCO (square; Refs. § 
andllq). The solid (red) circles denote the same Tl-2201 sam- 
ple measured at two different doping levels, as discussed in the 
text. Spectroscopic measurements of the gap from ARPES on 
Tl-2201 (red down-triangle; Refs. E3 and!) and YBCO (blue 
down triangle; Refjjjg) and SIS tunneling on Bi-2212 (black 
up-triangles; Ref. [191 ) are also shown. The dashed line gives 
the BCS doping dependence of the gap based on the mea- 
sured T c , namely Ao = uksTc, using a strong-coupling value 
of fi = 3 and T c (p) from Eq. [T] Lower panel: The gap-to- T c 
ratio, Ao/fcsTc, vs p, for Tl-2201 from thermal conductivity 
data. The solid (red) line is a fit through the two data points 
(full red circles) obtained on a single sample, for which the 
relative uncertainty is minimal (~ 3 %), given the fixed geo- 
metric factor. The dashed (red) lines on either side reflect the 
uncertainty on the absolute value (±17 %; see Table I). The 
horizontal (black) dashed line is the weak-coupling d-wave 
BCS prediction, namely A /k B T c = 2.14M 



1 dA 



hkf 



(4) 



proportional to the slope of the gap at the node. Note 
that, unique amongst transport properties, kq/T was 
shown to be robust against both Fermi liquid and ver- 
tex corrections^ 



The above formalism, however, is only valid in the limit 
ksT <C 7 <C Ao- Before proceeding, it is important to 
determine whether this criterion is indeed satisfied for our 
samples, particularly since Ao decreases with increasing 
doping and eventually goes to zero at the superconduc- 
tor/metal phase transition. — In the Appendix, we pro- 
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vide two independent and consistent estimates of 7 for 
each sample which show that in all cases 'y/ksT > 10 
(for T < 1 K) and 7/A0 < 0.16, validating the use of 
Eq.EJfor our Tl-2201 samples. 

The significance of Eq. [3] is that it provides a reli- 
able way to measure the gap from a rather straightfor- 
ward measurement of the thermal conductivity. Indeed, 
through Eqs. [3] and [U the slope of the gap at the node, 
6A/5cf)\ no de, is immediately obtained from kq/T, pro- 
vided Vp and kp are known* 6 .^ Making the assumption 
that A has the standard d x 2_ y 2 form A(<fi) — A cos(20), 
Eq. [3] can be rewritten (for vp 3> v a) as 



k% n v F 
6 c F A ' 



(5) 



kq/T can thus be related in a straightforward fashion to 
the gap maximum, Ao, by making use of the fact that 
vp and kp at the nodes are approximately doping inde- 
pendent. ARPES measurements 2 ^ have shown vf to be 
independent of both doping and material and to have a 
value vf = 2.7 xlO 7 ± 20% cm/s. Measurements of the 
Fermi surface by ARPES in LSCO 23 , Bi-2212^ and Tl- 
2201—, show Uf to have a value ~ 0.7 A^ 1 along the 
nodal direction (0,0) - (77,71"), essentially independent of 
material or doping {Uf only has significant doping de- 
pendence along (0,0) - (0,7r)). In overdoped Tl-2201, this 
value of kp was nicely confirmed by bulk sensitive angular 
magneto-resistance oscillation (AMRO) measurements* 2 ^. 
The values of Ao obtained for our five Tl-2201 samples 
from Eq. [5] are listed in Table I, and plotted vs p in the 
top panel of Fig. [3j Also shown is the value of Ao for an 
ultra-pure overdoped YBCO crystal^^ obtained in the 
same way (with Ko/T = 0.16 mW K~ 2 cm -1 ) - it is in 
excellent quantitative agreement with the Tl-2201 data. 

Ao decreases with increasing doping and follows the 
doping dependence of the gap measured by spectroscopic 
probes such as tunneling 19,26 and ARPES 27 . In Figure [3J 
we also plot the doping dependence of Ao in Bi-2212 mea- 
sured using SIS tunneling^ and find good quantitative 
agreement with Ao measured by thermal conductivity^ 
similar to previous work on optimally doped Bi-2212. 6 
Comparing directly to measurements on Tl-2201, re- 
cent ARPES measurements give Ao — 12 ± 5 meV and 
A ~ 24 ± 9 meV in a samples with T c = 30 K and 74 K 
respectively 9 ^ 7 - which is also in favorable agreement with 
thermal conductivity. (Ao here is taken as the average 
of two commonly used measures of Ao; the peak position 
and the leading-edge midpoint, which are 17 meV (33 
meV) and 8 meV (15 meV), respectively, for T c = 30 K 
(74 K).) 

The agreement between the different measures of A has 
important implications. First, it provides further com- 
pelling evidence for coherent nodal quasiparticle excita- 
tions that are correctly described by mean-field rf-wave 
quasiparticle theory throughout the overdoped region of 
the phase diagram. Secondly, the fact that the gap de- 
rived from thermal conductivity (i.e. from low-energy 



nodal quasiparticles) agrees well with the gap maximum 
derived from tunneling or ARPES (determined from anti- 
nodal quasiparticles at higher energies) suggests that the 
gap structure (as a function of angle) is not very differ- 
ent from the simple d-wave gap assumed in our analysis 
(Aocos(20)). This places limits on any sub-dominant 
order parameter that may exist in the bulk. Subdomi- 
nant order parameters that result in a fully gapped nodal 
region (such as is or id) are limited to having a magni- 
tude much less than the impurity bandwidth 7^ which 
is estimated to be on the order of Ao/10 in our Tl-2201 
samples (see Appendix) . An even more stringent restric- 
tion is provided by the YBCO measurements, where, in 
crystals of the highest quality, 7 can be orders of magni- 
tude smaller than in Tl-22014^ These conclusions are 
in agreement with the phase-sensitive tricrystal experi- 
ments, which show d x 2_ y 2 symmetry of the order param- 
eter and no signature of a subdominant imaginary con- 
tribution to the superconducting order parameter away 
from optimal doping. 30 (Note that a real s component 31 
is allowed.) 

The gap-to-T c ratio, fi = Ao/fcsTc, is plotted vs p in 
Figure [3J /x is seen to decrease with increasing doping 
from a strong-coupling BCS value of \i = 5-6 near opti- 
mal doping towards the weak-coupling value of fi = 2.14 
at high doping. This trend is made compelling by com- 
paring measurements on the same sample at two different 
dopings (solid red circles), for which the relative uncer- 
tainty is small. Indeed for sample a (see Table I), fi 
drops by a factor 1.9 (from 5.6 to 3.0) when the dop- 
ing is increased from p = 0.19 (T c = 84 K) to p = 0.25 
(T c = 26 K). 

Having established how the superconducting gap grows 
with decreasing doping in the overdoped regime of the 
phase diagram, it is is natural to ask how it continues to 
evolve below optimal doping, into the enigmatic under- 
doped regime. Heat transport measurements on YBCO 
have shown that kq/T decreases monotonically, at least 
down to p ~ 0.1i & 32 i 33 Based on Eq. 5, this was inter- 
preted in terms of a gap A that keeps on rising as p 
is reduced^ Note, however, that recent measurements 
of the in-plane resistivity in the (field-induced) normal 
state of YBCO reveal a "metal-insulator" crossover at 
p ~ 0.09.— One should not expect the standard BCS 
analysis (leading to Eq. 3) to apply below that crossover. 
(The same situation was encountered in LSCO^ set- 
ting in at a considerably higher doping since the "metal- 
insulator" crossover in that system is at p ~ 0.17M-) The 
crossover to a non-metallic normal state may account for 
the qualitative change observed in the behaviour of kq /T 
when p is reduced to 5-6%, where it becomes indepen- 
dent of both doping and magnetic field. 36 This overall 
evolution of kq/T vs p in YBCO is corroborated by re- 
cent ARPES measurements on Bi-2212, which show the 
gap near the nodes to have comparable magnitude and 
to rise with decreasing p until p ~ 0.08, and then fall 
abruptly* 3 - 7 - 

In conclusion, we have measured the doping depen- 
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dence of the low-temperature thermal conductivity in 
the overdoped cuprate Tl-2201. Analyzing the magni- 
tude of «o /T using mean- field theory, we show quantita- 
tive agreement between the slope of the gap at the nodes 
from thermal conductivity and spectroscopic measure- 
ments of the gap maximum at the anti-nodes in over- 
doped cuprates. These measurements provide further 
confirmation of the applicability of mean-field SCTMA 
transport theory in the cuprates, and more generally the 
basic description of overdoped cuprates as d-wave super- 
conductors. As a function of reduced doping, the gap- 
to- T c ratio increases by nearly a factor 2 from p = 0.25 
to p = 0.19. This reveals that the strong suppression 
of phase coherence (or T c ) characteristic of the under- 
doped regime, presumably caused by a combination of 
pair-breaking and phase fluctuations, in fact begins in 
the overdoped region, at least as high as p ~ 0.2, per- 
haps higher. 



APPENDIX A: ESTIMATES OF THE IMPURITY 
BANDWIDTH 

In this Appendix we estimate the impurity bandwidth, 
7, using two approaches: from the normal state trans- 
port and from the magnetic field dependence of kq /T in 
the superconducting state. Both approaches yield similar 
values of 7 which in all samples satisfy the two criteria 
for the validity of a "universal" limit analysis, namely 
7 < A and T < 7. 



1. From normal state transport 

First, we estimate 7 from the normal state quasipar- 
ticle scattering rate, IV . Because the normal state of 
our samples could not be reached with the magnetic 
field available in our lab (15 T), we turn to a previous 
study 7 on a crystal with a similar quality but signifi- 
cantly lower T c (higher doping), namely (15 K). In that 
crystal, Tjv was estimated from kn/T to be 0.52 meV* 7 . 
Using this value as representative of Tn in our samples, 
we can estimate an upper bound on 7 via the relation 
7 ~ 0.63v / Ia'Ao, valid in the unitarity limit (which 
yields the largest 7 for a given value of Tn)£& These 
estimates are listed in Table UH based on the A values 
of Tabic |H One can see that -f/k B T > 13 for T < 1 K 
and 7/A0 ~ 0.1 in all cases, more precisely ranging from 
0.07 to 0.16. When 7 is not negligible, there is a known 
correction to Eq. [3] (and Eq. |5j), which for 7/A0 ~ 0.1 is 
about 20%^ i.e. the clean-limit gap value is 20% larger 
than estimated using Eq. [5) An error of this magnitude 
has little impact on the conclusions of this article. 



TABLE II: Estimates of the impurity bandwidth, 7, from the 
normal state scattering rate Fjv = 0.52 meV, in the unitary 
limit. 
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2. From the magnetic field dependence of kq/T 

As a crosscheck, we can also make use of the magnetic 
field dependence of the thermal conductivity in the T — > 
limit to estimate 7 in our samples. Unlike the zero- 
field thermal conductivity, which is independent of 7, the 
field dependence of the thermal conductivity is dependent 
upon the quasiparticle scattering ratei 40 i 41 In a d-wave 
superconductor kq/T increases with an applied magnetic 
field due to the Volovik effect 4 ^ and in optimally doped 
cuprate a 13 ! 41 ' 43 the thermal conductivity is well described 
by semi-classical calculations of the field dependence of 
kq/T based on the Volovik effect i 40 ' 44 




T[K) 

FIG. 4: (colour online) Thermal conductivity in an applied 
field (H//c) for the Tl-2201 sample with T c = 76 K, plotted 
as k/T vs T. ko/T increases by a factor of ~ 3 in an applied 
field of 11.5 T. The lines are fits to Eq.[2] Inset: The in-field 
electronic thermal conductivity Ko(H)/T normalized to the 
zero-field value ko(0)/T as a function of magnetic field in a 
Tl-2201 sample with T c — 89 K (not used in the main study). 
The solid line is a fit to the semi-classical model (Eq. lAlj l. 
which yields 7/A0 — 0.1 (see text). 

An analytic model of kq(H)/T for H <C H C 2 and J _L 
H by Kiibert and Hirschfeld 40 gives: 

^= —f . (Al) 

ko/ 2 py/l + p 2 - smb. p 
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where kq/T is the universal value in zero field (Eq. [3]) 
and (for 7 > E H ) 

p = ^%hi 1 /E H . (A2) 

Here Eh = ah-J2pKVp •v/iJ/^o is the average energy 
shift experienced by a quasiparticle due to the Volovik 
effect, a is a constant dependent upon the geometry of 
the vortex lattice (= 1/2 for a square lattice 4 ^) and $0 
is the quantum of flux. Note that aside from 7, Eq. I All 
is only dependent upon known constants. As such, we 
can estimate 7 from a single parameter fit to the field 
dependence of kq/T at H <C H c2 . 

In Figure [U we show k/T vs T for the Tl-2201 sample 
with T c = 76 K, in H = T and 11.5 T. We see that n /T 
is enhanced by a factor 3 or so. Similar field-induced en- 
hancements of kq/T are observed in all samples, by a 
factor which varies from 2 to 44£ In the inset, we show 
Ko(H)/T (normalized to kq(0)/T) vs H in a separate 
Tl-2201 sample (T c = 89 K). 47 Eq. ED provides a reason- 
able fit (solid line) to the data, giving p\fii — 3.36T 1 / 2 , 



which results in 7 = 3.8 meV, so that 7/A0 — 0.1. This 
is similar in magnitude to our above estimate from nor- 
mal state transport, providing a non-trivial validation for 
both methods of estimating 7. 

We conclude that within possible corrections on the 
order of 20% or so, it is legitimate to use Eqs. OandOto 
extract an estimate of Ao from a measurement of kq /T 
in the Tl-2201 samples considered here. 
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